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We develop a bosonization scheme for the collective dynamics of a spinless two-dimensional electron 
gas (2DEG) in the lowest Landau level. The system is treated as a continuous elastic medium, and 
quantum commutation relations are imposed between orthogonal components of the elastic displace- 
ment field. This theory provides a unified description of bulk and edge excitations of compressible 
and incompressible phases, and explains the results of recent tunneling experiments at the edge of 
the 2DEG. 
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Bosonization M is a powerful mathematical technique 
to describe the long-wavelength dynamics of correlated 
electronic systems, in which the motion of the electrons 
is highly collectivized. This technique has proved very 
useful in studies of the behavior of interacting electrons 
in one-dimensional metals. Higher-dimensional general- 
izations have also been proposed g] . More recently, it has 
been realized that the edges of a two-dimensional electron 
gas (2DEG) in a strong magnetic field behave as a "chi- 
ral Luttinger liquid" H and are therefore amenable to 
bosonization HQ. A collective description of the uni- 
form incompressible 2DEG in a strong magnetic field 
such that all the electrons reside within the lowest Lan- 
dau level (LLL) || was pioneered by Girvin et al. 
and an effective bosonic theory was developed shortly af- 
terwards 0. The similarity between the local structure 
of the 2DEG in the LLL and that of a classical Wigner 
crystal was exploited by Johansson and Kinaret || to 
set up an "independent boson model" for the electronic 
spectral function; an equivalent model was derived by 
Haussmann || from diagrammatic theory. Finally, sev- 
eral authors have used bosonization schemes to describe 
the dynamics of a smooth compressible edge in a quan- 
tum Hall liqui d p~0| — , or of the electron gas in a weak 
magnetic field |13|. 

In this Letter we formulate a bosonic theory of the lin- 
ear dynamics of a general uniform or nonuniform distri- 
bution of electrons in the LLL. Unlike previous theories 
|0-|l3| , the present one treats bulk and edge excitations of 
compressible and incompressible states on an equal foot- 
ing. One of our main results is a simple explanation for 
recent experimental observations of one-electron tunnel- 
ing characteristics in the edge Jl4|] of the 2DEG. 

The basic idea of our approach is to treat the electrons 
in the LLL as a continuous elastic medium characterized 
by an equilibrium density po(r), and by a local displace- 
ment field u{r,t), such that the time-dependent density 
is (to first order in u) p(f, t) = po(f) — V • [pa{r)u(f, t)\ = 
Poif) + 5p(f, t). We ignore the spin in this paper. Under 



the assumption that both the equilibrium density and 
the displacement are slowly varying on the scale of the 
magnetic length I — (hc/eB) 1 / 2 (where B is the magnetic 
field) the effective long-wavelength Hamiltonian is 

H=\ Jdrjdr'5p{r) I - x~ X (f> ?) J 5p{?) 



dr/x \p (r*)] 

a, (3 



Sal3{r) - -5 a f3V ■ u{r) 



(1) 



where x _1 (r, r 1 ) is the inverse of the proper static density- 
density response function flq| , p, is the shear modulus 
(discussed below), s a p = [du a {f) / drp + du^r) / dr a }/2 is 
the strain tensor, and a, (5 are cartesian indices. Eq. (|l|) 
generalizes the Hamiltonian of classical elasticity theory 
]l6t to take into account the long range of the Coulomb 
interaction, and of x(r, r 1 ) in the incompressible states of 
the 2DEG §]. 

The second term of Eq. (0) is the crucial new feature 
of our theory. Ordinarily, the shear modulus of a liquid is 
taken to be zero, but this is only true for shear deforma- 
tions which occur slowly compared to the thermalization 
time t of the system. At frequencies higher than 1 /r the 
correlated electron liquid behaves as a Wigner glass, and 
therefore exhibits a nonvanishing shear modulus. The es- 
sential assumption is that, in this system, the relaxation 
frequency scale 1/r is much smaller than the character- 
istic frequency of the collective dynamics. 



The algebra of the displacement operators u x , 



'li- 



the LLL, is deduced from the canonical quantization con- 
dition for the hydrodynamical momentum and displace- 
ment fields \p a {r), up(r')] = —i6(r — r')6 al 3, by projecting 
out the higher Landau levels, i.e. averaging over the fast 
cyclotron motion. This leads to 



[u a (r),up(r')] = -ie a p5(r - r') I 2 / Po(r) , 



(2) 



where e a p is the two-dimensional Levi-Civita tensor. 
This is consistent with the commutation relation between 



1 



projected density fluctuations ra| in the long wavelength 
limit, namely [p(f),p(r')] = — iV r x V r < po(r)l 2 5(r — r'). 
Bosonization of p) is accomplished by 



u a {r) 



I 



1/2,- £ [ b n9na{r) + ^5«a(r)] 
Po ( r ) n>0 



(3) 



where 6„ and fct are boson operators, satisfying the 
usual algebra [&„,&„,] = S nn > etc.. The functions g na (r) 
((7* Q (r*)), with n > 0, are the positive (negative) fre- 
quency solutions of the eigenvalue problem 



Ha/3{f,T J )g n p(r')dr' = uj n g na (r) 



(4) 



which is equivalent to the classical equation of motion. 
(We adopt the convention of summing over repeated 
cartesian indices.) The operator H a p is given by 

H a p(?,?) = -il 2 P y\r)e ai [d^x-H^r')] p^P) 

I 2 r - i 

+ij-^e a ~ f [dpd' 1 -d 1 d' l3 + 8 1 pd n d' 71 ] 5(f-r')p, , (5) 

where x~ (r,r') = x" 1 ^,?"') — e 2 / 1 r^* — r-' | , and d a and dp 
denote derivatives with respect to r a and r'g respectively. 
It is hermitian with respect to the scalar product (/, g) = 
J /a(^)* e Q/3fl pffidr, and its eigenvalues are therefore real. 
For each eigenfunction g n with positive eigenvalue u> n , 
there is a complex conjugate one g* = g_ n with negative 
eigenvalue — uj n . These eigenf unctions can be chosen to 
satisfy orthonormality relations (g n ,9m) — sgn(n)S nm (n 
and m can now have either sign) JlTj . The "completeness 
relation" has the form 

^2\9 na (r)9np(r') - 9na(r)9n0(r')] = ie af) 8{r- ?). (6) 

n>0 

These equations guarantee that the commutation re- 
lation (0) is satisfied. Then, substitution of Eq. (||) 
into Eq. ([!]) yields the Hamiltonian in the desired form 

^ = En>o(^n + V2)w«. 

Our next step is to construct the electron tunneling 
operator ip(R) associated with the LLL coherent state 
orbital centered at R. This operator must satisfy the 
commutation relation 



mR),p(r)} = rp(R) exp[-|f - R\ 2 /2l 2 ]/2nl 2 



(7) 



i.e., destroy a gaussian density near R. Eq. (|7|) alone 
does not uniquely determine ip(R). In order to do this 
we must also specify its commutator with the vorticity 
p v (r) = — £ a (3d a [po{r)u[3(f)]. For the purpose of study- 
ing tunneling, we choose [ip(R), p v (r)] = 0, because the 
incoming electron is not expected to change the vorticity 
at the initial time. The solution has the form 



tp(R) = exp 



n>0 



M*(R) h M n (R) 



-b„ 



bt 



(8) 



where the "electron-phonon" matrix elements M n (R) are 
written as M n (R) = J exp[-|f- R\ 2 /2l 2 ]M n {r)df/2nl 2 , 
and 



M n (r) 



2^ I €a0 



pl /2 {?)l 



dp log 



\dr' 



(9) 



The calculation of the local Green's function G(R, t) = 
—i(Tip(R,t)ilj ji (R,0)} can be carried out by standard 
techniques JT^,^. Here we simply report the zero tem- 
perature result for the integral equation connecting the 
electronic spectral function A(R,lu) (the Fourier trans- 
form of G(R,t)) to the collective excitation spectrum: 



(10) 



ujA(R,lj)= g(R,ri)A(R,u; - fl)dQ 
Jo 

where uj > and 

g(R, ^) = ^ £ \M n (R)\ 2 S(n - uj n ) 



(11) 



n>0 



is closely related to the local dynamical structure factor 
of the liquid. The above equations constitute a complete 
scheme for the calculation of single particle and collective 
properties of a general distribution of electrons in the 
LLL. We now discuss some specific examples. 

(1) The uniform 2DEG. Let po be the uniform density. 
The normalized eigenfunctions of Eq. (|]) have the form 



9qL{r) = iql 



2uj q po 



1/2 



and 



where 



r -« 1 ( u q Po 



1/2 



P 



1/2 



v(q) 



K(q)+p 



Pi 



1/2 



(12) 



(13) 



(14) 



is the frequency of the mode at wavevector q, v(q) is 
the Fourier transform of the electron-electron interaction, 



K{q) 



-PoX 



_1 (g) is the g-dependent bulk modulus, and 



the labels L (longitudinal) and T (transverse) refer to 
the components parallel and perpendicular to q. The 
electron-phonon matrix element of Eq. (pt) has the form 



M q (r) 



v{q) 



KM 



P 



Pi 



p q (?) , (15) 



where p q = -iqlp)j 2 g qL (r) and K?{q) = -p%J x l (r- 
r')exp[z<f ■ (r' — r)]dr' . (The reason for the apparently 
unnecessary subscript r will become clear below). 

We now distinguish two cases: (i) Compressible case. 
limg^o K{q) = K is finite. The values of K and p can 



2 



be approximated by those of a classical Wigner crystal, 
namely, fj, ~ [i c = 0.09775p (e 2 /l)is 1/2 and K ~ K c = 
6/i c [jljj. The long wavelength modes have frequencies 
u} q oc g 3 / 2 for Coulomb- like interaction (v(q) = 2ne 2 /q) 



and 



q for short-range interaction (v(q) 



>d). 



We have now all the necessary input for calculating the 
spectral function from Eqs. ([l^jl^). The low- frequency 
behavior (where "zero" frequency is the chemical poten- 



tial) is found to be A(u>) oc uj 5 / 4 e 7 / w 1 for Coulomb 
interaction and A(w) oc w <?+W) 1 "*/l) 1 "/4>'-i for short _ 
range interaction. Complete numerical results are shown 
m Fig. @. They are in good agreement with previous cal- 
culations as well as with experimental data [pl| . 

(it) Incompressible case. Within the single-pole ap- 
proximation K(q — » 0) ~ Apo/2ag 4 , where A is the 
collective excitation gap at q = 0, and a = (1 — v)j%v 
is derived from the Laughlin wavefunction at filling fac- 
tors v = 2tt1 2 pq = 1/odd integer ||. Using Eq. (|l4| ) we 
find at long wavelength uj q = (A/i/2apo) 1 ^ 2 independent 
of q. This formula can be used to deduce the value of 
\x if A = lim 9 _»o is known. Alternatively, one can 
substitute \i — [i c from the classical Wigner crystal and 
obtain A„ = 0.391^ 3 / 2 (1 - i/pV/Z, which gives A 1/3 = 
O.lle 2 //, A 1/5 = 0.044e 2 // and A 1/7 = 0.025e 2 //. These 
results are in good agreement with variational estimates 
H^2|, with the exception of the first one which is al- 
most 30% lower than the variational one, but compares 
favorably with exact diagonalization studies Q . 

The calculation of the spectral function is more sub- 
tle. Straightforward application of the linear response 
formulas ( |10|{15| ) is incorrect, because the addition (or re- 
moval) of charge at point R creates a compressible region 
in the middle of the liquid, changing the topology of the 
incompressible region from simply to doubly connected. 
The change in topology can be taken into account in the 
following manner. We stipulate that the bulk modulus 
kernel — /OqX _1 (^, f') has the form characteristic of an in- 
compressible liquid (namely, \r — r'| 2 log|r — r'\) when 
both r and r' are within the incompressible region, but 
it is given by the local form K c 5{r — r') when either r 
or r' arc within the compressible "core" of the excita- 
tion. Because the size of the core region is microscopic, 
its presence does not affect significantly the frequencies 
of the long wavelength modes. On the other hand, in the 
"electron-phonon" matrix element, given by Eq. (|l^), we 
must use Kp{q) ~ K Cl because r is inside the core re- 
gion. Within this scheme, the calculation of the spectral 
function can be straightforwardly carried out, without 
adjustable parameters. In Fig. |l|, we plot the results for 
v = 1/3. We have used ^-dependent K and /i in or- 
der to fit accurately the collective mode dispersion and 
structure factor || at finite q. The essential difference be- 
tween this and the compressible case is the appearance 
of a 5- function peak at u) = 0, which now corresponds 
to fi + (recall that in the incompressible liquid the chem- 
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FIG. 1. Spectral function for a compressible liquid at 
v — 0.3 with short range interactions (d — 101, dashed curve), 
with long-range interactions (dotted curve) and for an incom- 
pressible liquid at v = 1/3 (full curve, the 5- function peak at 
u — contains around 2% of the spectral strength) . 



ical potentials for addition or removal of charge, /i+ and 
are different pi)] ). The strength of the 5-function 
is Z = exp[— J2 q \M q 1 2 /u>g] , which would have vanished 
in the compressible case, and does not vanish here be- 
cause of the gap. The peak at /x+ reflects the ability of 
the incompressible liquid to accommodate the incoming 
electron in the ground state as topological defects of the 
initial incompressible state, without creating collective 
excitations. The incoherent part of the spectral func- 
tion, at higher frequencies, corresponds to the creation 
of additional collective excitations |25) . 

(2) Edge dynamics. The case of a smooth compressible 
edge has been treated in previous publications |ll| . Here 
we focus on the case of a sharp edge, which is directly 
relevant to the interpretation of recent lateral tunneling 
experiments by Chang et al. and Grayson et al. [[Tif . Let 
us consider a straight edge along the y axis, and let eEx 
be the gradient of the confinement potential at the edge. 
The density po is assumed to be uniform for x < 0, and 
zero for x > 0. The presence of the edge electric field 
breaks the rotational symmetry of the Hamiltonian (Q) , 
and must be taken into account with the additional term 



^jedge 



e / p(r)u(r) ■ E dr . 



(16) 



With this term included into the eigenvalue problem (0), 
we obtain a new set of solutions, which satisfy the condi- 
tions V-g = and Vxg = for x < 0. Neglecting the log- 
arithmic correction arising from the long-range Coulomb 
interaction, they obey the equation of motion 



iug a (r) = ve a pd{ig x (f) , 



(17) 
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where v = cE/B is the classical drift velocity. The or- 
thonormal solutions are (for x < 0) 



9q(r) 



= q 1 / 2 e qx - iqv {x 



iy) 



(18) 



where q > is a one-dimensional vector along the edge, 
and the eigenvalues are uj q = vq. These are analogous to 
gravity waves on the surface of a liquid. Because there is 
neither density change nor shear strain in the interior of 
the system, these solutions do not depend on the values 
of the bulk elastic constants. The use of elasticity theory 
is justified at small q since the displacement field u is 
slowly varying. 

The effective edge dynamics can be derived from the 
full dynamics by projecting the latter onto the subspace 
spanned by the edge- wave solutions of Eq. (|l8|). Within 
this subspace, we can define an "edge density" operator 
p edge (y) = J Sp(x,y)dx = pou x (0,y), where the bar de- 
fer 



notes projection onto the edge-wave subspace, i.e. 



instance, u x (0,y) = [l/p]/ 2 ) E„> <? 1/2 [& 9 e 



g>0 



It is easy to verify that the edge density satisfies the stan- 
dard Kac-Moody algebra § [pf gc ,p c ^\ = (vq/2ic)6 qq , 
where v is the bulk filling factor. Thus, we have deduced 
the dynamics of the chiral Luttinger liquid at the edge 
from a projection of the canonical dynamics of displace- 
ment fields in the bulk. The edge tunneling operator, 
obtained by imposing [tp(y), p edse (y')] — ip(y)S(y — y') in 
the edge- wave subspace, is still given by Eq. (||), with 
the sum running over the edge modes, and the matrix el- 
ement M q {y) = v{2TTq/v) 1 ^ 2 e- iqv . Use of Eqs. (jTo)) and 
( pi] ) then leads to the conclusion that the tunneling cur- 
rent must vanish, at low bias V, as V 1 ^ . This result is 
in good agreement with the experimental findings. The 
present derivation explains why the tunneling exponent 
is found to depend only on the bulk density, and not on 
whether the bulk is compressible or not. 

In summary, we have developed a magneto-elastic 
bosonization scheme for the long wavelength dynamics of 
the 2DEG in the LLL. Our results show that this scheme 
can provide a unified description of different physical ef- 
fects in the bulk and at the edge of the system. 

This work was supported by NSF grant No. DMR- 
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L. Glazman for a useful discussion at the beginning of 
this work. 
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